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Abstract 

We rigorously prove the existence of chaotic dynamics for the triopoly game 
model already studied, mainly from a numerical viewpoint, in [To]. In the model 
considered, the three firms are heterogeneous and in fact each of them adopts 
a different decisional mechanism, i.e., linear approximation, best response and 
gradient mechanisms, respectively. 

The method we employ is the so-called "Stretching Along the Paths" (SAP) 
technique in [20J, based on the Poincare-Miranda Theorem and on the proper- 
ties of the cutting surfaces. 

Keywords: Chaotic dynamics; Stretching along the paths; triopoly games; hetero- 
geneous players. 

2000 AMS subject classification: 54H20, 54H25, 37B10, 37N40, 91B55. 

1 Introduction 



In the economic literature, due to the complexity of the models considered, an 
analytical study of the associated dynamical features turns out often to be too difficult 
or simply impossible to perform. That is why many dynamical systems are studied 
mainly from a numerical viewpoint (see, for instance, [21 IH EDI EI])- Sometimes, 
however, even such kind of study turns out to be problematic, especially with high 
dimensional systems, where several variables are involved. 

In particular, as observed in Naimzada and Tramontana's working paper [15], this 
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may be the reason for the relatively low number of works on triopoly games (see, for 
instance, [9j [23j [29]), where the context is given by an oligopoly composed by three 
firms. In such framework, a local analysis can generally be performed in the special 
case of homogeneous triopoly models, i.e., those in which the equations describing 
the dynamics are symmetric (see, for instance, [T| 13"! I2"l]). 

A more difficult task is that of studying heterogeneous triopolies, where instead the 
three firms considered behave according to different strategies. This has been done, 
for instance, in [11], as well as in the above mentioned paper by Naimzada and 
Tramontana [TH] where, in addition to the classical heterogeneity with interacting 
agents adopting gradient and best response mechanisms, it is assumed that one of 
the firms adopts a linear approximation mechanism, which means that the firm does 
not know the shape of the demand function and thus builds a conjectured demand 
function through the local knowledge of the true demand function. In regard to such 
model, those authors perform a stability analysis of the Nash equilibrium and show 
numerically that, according to the choice of the parameter values, it undergoes a flip 
bifurcation or a Neimark-Sacker bifurcation leading to chaos. 

What we then aim to do in the present paper is complementing that analysis, by 
proving the existence of chaotic sets only via topological arguments. This task will 
be performed using the "Stretching Along the Paths" (from now on, SAP) technique, 
already employed in [T3] to rigorously prove the presence of chaos for some discrete- 
time one- and bidimensional economic models of the classes of overlapping generations 
and duopoly game models. Notice however that, to the best of our knowledge, this is 
the first three-dimensional discrete-time application of the SAP technique, called in 
this way because it concerns maps that expand the arcs along one direction. We stress 
that, differently from other methods for the search of fixed points and the detection 
of chaotic dynamics based on more sophisticated algebraic or geometric tools, such 
as the Conley index or the Lefschetz number (see, for instance, [El Q31 128] ) , the SAP 
method relies on relatively elementary arguments and it is easy to apply in practical 
contexts, without the need of ad-hoc constructions. No differentiability conditions 
are required for the map describing the dynamical system under analysis and even 
continuity is needed only on particular subsets of its domain. Moreover, the SAP 
technique can be used to rigorously prove the presence of chaos also for continuous- 
time dynamical systems. In fact, in such framework it suffices to apply the results in 
Section [21 suitably modified, to the Poincare map associated to the considered system 
and thus one is led back to work with a discrete-time dynamical system. However, 
the geometry required to apply the SAP method turns out to be quite different in the 
two contexts: in the case of discrete-time dynamical systems we look for "topological 
horseshoes" (see, for instance, [3 [121 [32] ) , that is, a weaker version of the celebrated 
Smale horseshoe in [27J, while in the case of continuous-time dynamical systems one 
has to consider the case of switching systems and the needed geometry is usually that 
of the so-called "Linked Twist Maps" (LTMs) (see [H1E1I22]), as shown for the planar 
case in [T51 l2"Tj . We also stress that the Poincare map is a homeomorphism onto its 
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image, while in the discrete-time framework the function describing the considered 
dynamical system need not be one-to-one, like in our example in Section |3j Hence, 
in the latter context, it is in general not be possible to apply the results for the 
Smale horseshoe, where one deals with homeomorphisms or diffeomorphisms. As 
regards three-dimensional continuous-time applications of the SAP method, those 
have recently been performed in [25J, in a higher-dimensional counterpart of the 
LTMs framework, and in (26], where a system switching between different regimes is 
considered. 

For the reader's convenience, we are going to recall in Section [2] what are the basic 
mathematical ingredients behind the SAP method, as well as the main conclusions it 
allows to draw about the chaotic features of the model under analysis. It will then 
be shown in Section [3] how it can be applied to the triopoly game model taken from 
[15] . Some further considerations and comments can be found in Section HJ which 
concludes the paper. 



2 The "Stretching along the paths" method 



In this section we briefly recall what the "Stretching along the paths" (SAP) 
technique consists in, referring the reader interested in further mathematical details 
to [20], where the original planar theory by Papini and Zanolin in [TBI [T7] has been 
extended to the N— dimensional setting, with N > 2. 

In the bidimensional setting, elementary theorems from plane topology suffice, while 
in the higher-dimensional framework some results from degree theory are needed, 
leading to the study of the so-called "cutting surfaces". In fact, the proofs of the 
main results in [20J (and in particular of Theorem 12. II below), we do not recall here, 
are based on the properties of the cutting surfaces and on the Poincare-Miranda 
Theorem, that is, an N- dimensional version of the Intermediate Value Theorem. 
Since in Section [3] we will deal with the three-dimensional setting only, we directly 
present the theoretical results in the special case in which N = 3. 

We start with some basic definitions. 
A path in a metric space X is a continuous map 7 : [to^i] X. We also set 7 := 
7([to,ti]). Without loss of generality, we usually take the unit interval [0,1] as the 
domain of 7. A sub-path a of 7 is the restriction of 7 to a compact sub-interval of its 
domain. By a generalized parallelepiped we mean a set V C X which is homeomorphic 
to the unit cube I 3 := [0, l] 3 , through a homeomorphism h : M. 3 D I 3 — y V C X. We 
also set 

VJ := h([x 3 = 0}) , V- := h([x 3 = 1]) 
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and call them the left and the right faces of V, respectively, wher^j] 

[x 3 = 0] := {(x 1 ,x 2 ,x 3 ) el 3 : x 3 = 0} and [x 3 = 1] := {(x 1 ,x 2 ,x 3 ) e I 3 : x 3 = 1}. 

Setting 

v~ ■■= vj u v~ , 

we call the pair 

V ■= (V,V~) 

an oriented parallelepiped ofX. 

Although in the application discussed in the present paper the space X is simply 
M 3 and the generalized parallelepipeds are standard parallelepipeds, the generality of 
our definitions makes them applicable in different contexts (see Figure 1). 

We are now ready to introduce the stretching along the paths property for maps 
between oriented rectangles. 

Definition 2.1 (SAP) Let A := {A, A") and B := (B,B~) be oriented paral- 
lelepipeds of a metric space X. Let also i|i : I be a function and /C C A 
be a compact set. We say that (/C, i/j) stretches A to B along the paths, and write 

(JC,ip) : A^B, (2.1) 

if the following conditions hold: 

• ij} is continuous on K, ; 

• for every path 7 : [0, 1] — > A with 7(0) and 7(1) belonging to different com- 
ponents of A~, there exists a sub-path a := 7 1 [?,*"] : [0, 1] D [t',t"] — > /C, such 
that i>{o-{t)) e B,Vt e [t',t"], and, moreover, 4>{cr(t')) and if)(a(t")) belong to 
different components of B~ . 

For a description of the relationship between the SAP relation and other "covering 
relations" in the literature on expansive-contractive maps, we refer the interested 
reader to [T3] . 

A first crucial feature of the SAP relation is that, when it is satisfied with A = B%, 
it ensures the existence of a fixed point localized in the compact set /C. In fact the 
following result does hold true. 

1 Notice that the choice of privileging the third coordinate is purely conventional. In fact, any 
other choice would give the same results, as it is possible to compose the homeomorphism h with a 
suitable permutation on three elements, without modifying its image set. 

2 Note that this means both that A and B coincide as subsets of X and that they have the same 
orientation. In fact, it is easy to find counterexamples to Theorem 12.11 if the latter property is 
violated (see, for instance, [T5], pag. 11). 
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Theorem 2.1 Let V := {V,V ) be an oriented parallelepiped of a metric space X 
and let ip : V — > X be a function. If K CP is a compact set such that 

(/C,V) :V*»V, 
then there exists at least a point z £ K, with ifi(z) = z. 

For a proof, see [20J, pagg. 307-308. Notice that the arguments employed therein 
are different from the ones used to prove the same result in the planar context (see, 
for instance, [13], pagg. 3301-3302), which are in fact much more elementary. 
A graphical illustration of Theorem 12.11 can be found in Figure 1, where it looks 
evident that, differently from the classical Rothe and Brouwer Theorems, we do not 
require that ip(dA) C A. 




Figure 1: The tubular sets A and B in the picture are two generalized parallelepipeds, 
for which we have put in evidence the compact set K. and the boundary sets Aj and 
A~ , as well as and B~ . In this case (JC, ip) stretches the paths of A across A itself 
and therefore the existence of a fixed point for %p in /C is ensured by Theorem 12.11 

The most interesting case in view of detecting chaotic dynamics is when there 
exist pairwise disjoint compact sets playing the role of K, in Definition 12.11 Indeed, 
applying Theorem 12.11 with respect to each of them, we get a multiplicity of fixed 
points localized in those compact sets. Another crucial property of the SAP relation 
is that it is preserved under composition of maps, and thus, when dealing with the 
iterates of the function under consideration, it allows to detect the presence of peri- 
odic points of any period (see Lemma A.l, Theorems A.l and A. 2 in [13] . which can 
be directly transposed to the three-dimensional setting, with the same proof). 
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We now describe in Definition 12.21 what we mean when we talk about "chaos" and 
we explain in Theorem 12.21 which is the relationship between that concept and the 
stretching relation in Definition 12. II We stress that Theorem 12.21 is the main theoret- 
ical result we are going to apply in Section [3] and that it can be shown exploiting the 
two properties of the SAP relation mentioned above. In fact, its proof follows by the 
same arguments in Theorems 2.2 and 2.3 in |13j . 



Definition 2.2 Let X be a metric space and let if) : X D D — » X be a function. 
Let also m > 2 be an integer and let /Co, . . . , /C m _i be nonempty pairwise disjoint 
compact subsets of V. We say that if) induces chaotic dynamics on m symbols on the 
set D relatively to /C , . . . , K, m -\ if, setting 



m—l 



(2.2) 



£ := [J £j C V 
and defining the nonempty compact set 

oo 
n=0 

then there exists a nonempty compact set 

X Q Zoo Q fc, 
on which the following conditions are fulfilled: 

(i) if}(X)=X; 

(ii) ip\x is semi-conjugate to the Bernoulli shift on m symbols, that is, there exists a 
continuous map tv : X — > S+, where := {0, 1} N is endowed with the distance 



m 



i+l 



for s' = (s-)ieN , s 



[Si JieN t ^„ 



is the discrete distance on {0,1}, i.e., d(s' i ,s") = for = s'/ and 
d( s 'n s i) — 1 f° r s 'i 7^ s i'), suc h that the diagram 



X 



X 



(2.3) 



commutes, where a : —> is the Bernoulli shift defined by cr((sj 
(s i+1 ) h Vz G N; 
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(Hi) the set of the periodic points of i\)\t x is dense in X and the pre-image vr _1 (s) C X 
of every fc-periodic sequence s = (sj)^ ^ ^1 contains at least one fc-periodic 
point. 

Remark 2.1 According to Theorem 2.2 in IT^ . from (ii) in Definition \2.2\ it follows 
that: 

— h top (ip) > ht op (i/;\z) > h top (a) = log(m), where h top is the topological entropy; 

— there exists a compact invariant set A CI such that is semi-conjugate to 
the Bernoulli shift on m symbols, topologically transitive and displays sensitive 
dependence on initial conditions. 

Theorem 2.2 Let V := (V,V~) be an oriented parallelepiped of a metric space X 
and let ip : V — ^ X be a function. If JCq, . . . , K, m -\ are m > 2 pairwise disjoint 
compact subsets of V such that 

(/Q,V) :V^V, /ori = 0,...,m-l, (2.4) 
then if) induces chaotic dynamics on m symbols on V relatively to JCq, ... , K, m -\- 

Notice that if the function if> in the above statement is also one-to-one on /C : = 
V)T=q then it is additionally possible to prove that ip restricted to a suitable 
invariant subset of /C is semi-conjugate to the two-sided Bernoulli shift a : S 2 — > S 2 , 
a((si)i) := (s i+1 )i, Vi E Z, where S 2 := {0, 1} Z (see [2H Lemma 3.2])|. 

We are now in position to explain what the SAP method consists in. Given a 
dynamical system generated by a map ip, our technique consists in finding a subset V 
of the domain of ip homeomorphic to the unit cube and at least two disjoint compact 
subsets of V for which the stretching property in (j2.4p is satisfied (when V is suitably 
oriented). In this way, Theorem 12.21 ensures the existence of chaotic dynamics in 
the sense of Definition 12.21 for the system under consideration and, in particular, the 
positivity of the topological entropy for ip, which is in fact generally considered as 
one of the trademark features of chaos. 

3 The triopoly game model 

In this section we apply the SAP method to an economic model belonging to the 
class of triopoly games, taken from [IB] . 

By oligopoly, economists denote a market form characterized by the presence of a 
small number of firms. Triopoly is a special case of oligopoly where the firms are 

3 This is not the case in our application in Section [H Indeed, as it looks clear from Figure 2, the 
map F in (|3.2[) is not injective on the set JCq U 1C\ introduced in Theorem I3.ll 
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three. The term game refers to the fact that the players - in our case the firms - make 
their decisions reacting to each other actual or expected moves, following a suitable 
strategy. In particular, we will deal with a dynamic game where moves are repeated 
in time, at discrete, uniform intervals. 

More precisely, the model analyzed can be described as follows. 

The economy consists of three firms producing an identical commodity at a constant 
unit cost, not necessarily equal for the three firms. The commodity is sold in a single 
market at a price which depends on total output through a given inverse demand 
function, known to one firm (say, Firm 2) globally and to another firm (say, Firm 
1) locally. In fact, Firm 1 linearly approximates the demand function around the 
latest realized pair of quantity and market price. Finally, Firm 3 does not know 
anything about the demand function and adopts a myopic adjustment mechanism, 
i.e., it increases or decreases its output according to the sign of the marginal profit 
from the last period. The goal of each firm is the maximization of profits, i.e., the 
difference between revenue and costs. The problem of each firm is to decide at the 
beginning of every time period t how much to produce in the same period on the basis 
of the limited information available and, in particular, on the expectations about its 
competitors' future decisions. 

In what follows, we introduce the needed notation and the postulated assumptions: 



1. Notation 

x t : output of Firm 1 at time t ; 
y t : output of Firm 2 at time t ; 
z t : output of Firm 3 at time t ; 
p: unit price of the single commodity . 



2. Inverse demand function 



p:= . (3.1) 

x + y + z 



3. Technology 



The unit cost of production for firm % is equal to q, i = 1, 2, 3, where c\, C2, C3 are 
(possibly different) positive constants. 



4. Price approximation 

Firm 1 observes the current market price p t and the corresponding total supplied 
quantity Q t = x t + y t + z t . By using market experiments, that player obtains the slope 
of the demand function at the point (Qt,Pt) and, in the absence of other information, 
it conjectures that the demand function, which has to pass through that point, is 
linear. 



8 



5. Expectations 

In the presence of incomplete information concerning their competitors' future 
decisions (and therefore about future prices), Firms 1 and 2 are assumed to use naive 
expectations. This means that at each time t both Firm 1 and 2 expect that the 
other two firms will keep output unchanged w.r.t. the previous period. 



As shown in [T5], the assumptions above lead to the following system of three differ- 
ence equations in the variables x, y and z: 

_ _ 2x t +yt+zt-c 1 (x t +yt+zt) 2 
•I't+l — 2 



Vt+l = yj ^ -Xt-Zt (TG) 

zt+i =* + azt (-eg + (x X+*ty< 

where a is a positive parameter denoting the speed of Firm 3's adjustment to changes 
in profit and c±, c 2 , c 3 are the marginal costs. 

We refer the interested reader to [15] for a more detailed explanation of the model, 
as well as for the derivation of (TG). 

As mentioned in the Introduction, in [15] Naimzada and Tramontana discuss the 
equilibrium solution of system (TG) along with its stability and provide numerical 
evidence of the presence of chaotic dynamics. In particular, it is shown the existence 
of a double route to chaos: according to the parameter values, the Nash equilibrium 
can undergo a flip bifurcation or a Neimark-S acker bifurcation. Moreover, in [T5] the 
authors numerically find multistability of different coexisting attractors and identify 
their basins of attraction through a global analysis. 

Hereinafter we will integrate that study rigorously proving that, for certain parameter 
configurations, system (TG) exhibits chaotic behavior in the precise sense discussed 
in Section 00. 

In order to apply the SAP method to analyze system (TG), it is expedient to represent 
it in the form of a continuous map F = (F l ,F 2 , F 3 ) : M. 3 + -> M 3 , with components 

F x {x,y,z) := gg±fe|fe±^ ; 



F 2 {x,y,z):=yJ^-x-z, (3.2) 
F 3 (x, y, z) := z + az \-c 3 A '^ L - 



(x+y+z) 2 J ' 

We prove that the SAP property for the map F is satisfied when choosing a 
generalized rectangle in the family of parallelepipeds of the first quadrant described 



4 Notice that, as we shall stress in Section HJ we only prove existence of an invariant, chaotic set, 
not its attractiveness. 
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analytically by 

TZ = Ufa, yi, Zi) := {(x,y,z) G M 3 : x e < x < x r , y £ < y < y r , z t < z < z r }, (3.3) 

with xi < x r , yt < y r , z e < z r and Xi, y^, z i7 i e {£, r}, satisfying the conditions in 
Theorem 13.11 

The parallelepiped TZ can be oriented by setting 

TZJ := [x e ,x r ] x [yt,y r ] x {z £ } and Tl~ := [x e ,x r ] x [y e ,y r ] x {z r } ■ (3.4) 

Consistently with [15], we choose the marginal costs as c\ = 0.4, C2 = 0.55 and 
C3 = 0.6. On the other hand, in order to easily apply the SAP method we need the 
parameter a to be close to 17, while in [T5] the presence of chaos is numerically proven 
for a around 8 0. The implications of this discrepancy will be discussed in Section HJ 




Figure 2: A possible choice of the parallelepiped TZ for system (TG), according to 
conditions (H1)-(H5). It has been oriented by taking as [ • ]~-set the union of the two 
horizontal faces TZJ and TZ~ defined in ( I3.4|) . In addition to F(TZj) and F(TZ~), we 
also represent the image set of two vertical faces of TZ. Notice that we used the same 
color to depict a set and its F-image set. 

Our result on system (TG) can be stated as follows: 

Theorem 3.1 If the parameters of the map F defined in (13. 2p assume the following 
values 

Cl = 0.4, c 2 = 0.55, c 3 = 0.6, a = 17, (3.5) 

5 As explained below, it would be possible to apply our technique with a lower value for a, at the 
cost of changing the parameter conditions in Theorem 13.11 and of making the computations in the 
proof much more complicated. However, it seems not possible to apply the SAP method to the first 
iterate of F when a is close to 8, which is the largest value considered in |15j . 
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Figure 3: This picture complements the previous one, by showing how the two vertical 
faces of TZ not considered in Figure 1 are transformed by the map F. Again, the same 
color is used to depict a set and its F- image set. 

then, for any parallelepiped TZ = TZ{x^, y iy zi) belonging to the family described in (13.31) . 
with Xi, yi, Zi, i € {£,r}, satisfying the conditions: 

(HI) z e = 0; 

(H2) x e + y e > z r > J^^fae + yi) - (a* + yi) > ; 

(H3) 2 ^-£-(x r + fr) ~ {X r + Vr)) > Z r ! 

{HA) ± - x r > y r + z r > - x £ > , ^ - x r > y e + z e , x r > ^ , 
2^ (1 - c x {yt + y r + z £ + z r )) >xi>0, \J^ L - {yt + zi) > x t ; 
{Hh) Xi + ze>-±, y r >^^-{x i + z e )>0, - {x r + z r ) > y t > , 

and oriented as in (\3.AL there exist two disjoint compact subsets /Co = Kq(JZ) and 
/Ci = JCi{TZ) oflZ such that 

(Ki, F) : TZ^TZ, for i = 0, 1. (3.6) 

Hence, the map F induces chaotic dynamics on two symbols on TZ relatively to /Co 
and K\ and displays all the properties listed in Theorem \2.B. 

Before proving Theorem 13 .![ we make some comments on the conditions in {HI)- 
{H5). First of all, notice that those conditions imply that xi+zg > and xi+ye+ze > 
and thus there are no issues with the definition of F on TZ %. We also remark that we 

6 Notice that, with our conditions on the parameters, it is immediate to check that also the 
functions we will introduce in the proof of Theorem 13 . 1 1 will be well defined, even when not explicitly 
remarked. 
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Figure 4: With reference to the parallelepiped TZ in Figure 2, reproduced here at a 
different scale, we show that the F-image set of an arbitrary path 7 joining in TZ the 
two components of the boundary set TZ~ intersects TZ twice. In particular, this is 
due to the fact that the horizontal faces TZJ and TZ~ are mapped by F below TZJ , in 
conformity with conditions (CI) and (C2), and that the flat surface S of the middle 
points w.r.t. the ^-coordinate in TZ is mapped by F above TZ~ , in agreement with 
condition (C3 ). 

chose to split (H1)-(H5) according to the corresponding conditions (C1)-(C5) in the 
next proof they allow to verify. Moreover we stress that the assumptions in {HI)- 
{H5) are consistent, i.e., there exist parameter configurations satisfying them all. For 
instance, we checked that they are fulfilled for C\ = 0.4, c 2 = 0.55, c 3 = 0.6, a = 17, 
x e = 0.5766666668, x r = 0.6316666668, y t = 0.3366666668, y r = .04516666668, z t = 
0, z r = 0.3951779684. These are the same parameter values we used to draw Figures 
2-6, with the only exception of zg that in those pictures is slightly negative. Although 
this makes no sense from an economic viewpoint, as the variables x, y and z represent 
the output of the three firms, we made such choice in order to make the pictures easier 
to read. In fact, choosing Z£ — 0, then F(TZj) C TZJ and thus the crucial set F(TZj) 
would have been not visible in Figures 2-4. With this respect, we also remark that 
in Figure 3 the x-axis has been reversed in order to make the double folding of F(TZ) 
more evident. 

In regard to the choice of the parameter values in (13.51) . as mentioned above, they are 
the same as in [15], except for a, which is larger here. In fact, numerical exercises 
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Figure 5: Since F(S) H TZ = (see Figure 4), then ft n F(K) =1C UK 1 , with /C 
and /Ci disjoint. 

we performed show that when a increases it becomes easier to find a domain where 
to apply the SAP technique. On the other hand, it seems not possible to apply 
our method for a sensibly smaller value of a. The impossibility of reducing a much 
below 17 comes from the fact that, as it is immediate to verify, when such parameter 
decreases it becomes more and more difficult to have all conditions in (H2) and (H3) 
fulfilled and with a = 10 it seems just impossible. The situation would slightly 
improve dealing with (C2) and (C3) below, instead of (C2) and (C3 ) as we actually 
do in order to simplify our argument, but still computer plots suggest it is not possible 
to have both conditions satisfied when a = 8, that is the largest value considered in 

mi- 

Proof. We show that, for the parameter values in ( I3.5p . any choice of Xi, Zi, i G 
{£,r}, fulfilling (H1)-(H5) guarantees that the image under the map F of any path 
7 = (71,72,72) : [0, 1] — > TZ = lZ(xi,yi, Zi) joining the sets TZJ and 7Z~ defined in 
( 13. 4 p satisfies the following conditions: 

(CI) F 3 ( 7 (0))<^; 

(C2) F 3 ( 7 (l))<^; 

(C3) 3 t*e (0,l):F 3 ( 7 (r))>^; 

(C4) ^(7^)) C [x t ,x r ], V* G [0,1]; 

(C5) F 2 ( 7 (t))C[ W) y r ],Vte[0,l]. 

Broadly speaking, conditions (C1)-(C3) describe an expansion with folding along the 
z-coordinate. In fact, the image F o 7 of any path 7 joining in TZ the sides TZJ and 
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Figure 6: Given the arbitrary path 7 in Figure 5 joining in TZ the two components 
of TZ~ , we show that the F-image sets of 7 fl /C and of 7 PI K,\ join TZj with TZ~ , as 
required by the SAP property 



TZ~ crosses a first time the parallelepiped TZ for t £ (0, £*) and then crosses TZ back 
again for t £ (£*, 1). Conditions ((74) and (C5) imply instead a contraction along the 
^-coordinate and the y-coordinate, respectively. 

Actually, in order to simplify the exposition, instead of ((73), we will check that the 
stronger condition 

(C3') F 3 (x,y,^) > z r , V{x,y) £ [x^x r ] x [y t ,y r ], 

is satisfied, which means that the inequality in ((73) holds for any t* £ (0, 1) such 
that j(t*) = (x, y, Xe 1 2 Xr ), for some (x, y) £ [xg, x r ] x [yi,y r ]- Notice that 

S := {(x,y, Ze + Zr ^ ; e [^,x r ] x [yi>,y r ]^ C TZ (3.7) 

is the flat surface of middle points w.r.t. the ^-coordinate in TZ depicted in Figure 4. 
Setting 

TZ := Ux,y,z) £ 1R 3 : (x,y) £ [ar/,ar r ] x \y e ,y r ], z £ 
:= Ux,y,z) £ R 3 : (x, y) £ [^,x r ] x [^,i/ r ], z £ 

and 

/C := TZ n F(^) and /d := fti fl F(ft) 

(see Figure 5), we claim that (CI), ((72), ((73'), ((74) and ((75) together imply (13TB]1 . 
Notice at first that /Co and /Ci are disjoint because, thanks to condition ((73 ), the 
set S in (13. 7p is mapped by F outside TZ (see Figure 4). Furthermore, by ((71), ((72) 



2* ? 



Z£ T" Z r 
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and (C3 ), for every path 7 : [0, 1] — > 1Z such that 7(0) and 7(1) belong to different 
components of 1Z~ , there exist two disjoint sub-intervals [to, to], [t'i,t'{] C [0,1] such 
that, setting a := 7l[^'] : [to, to] ~> ^0 and a x := j\[t' lt t'{] ■ ~> K-U it holds that 

F(cr (t )) and F(<7 (to)) belong to different components of 1Z~, as well as F(ai(t' l )) 
and F{ax{t'[)). Moreover, from (C4) and (C5) it follows that F(cr (t)) G K,Vt G 
[f ,tg] and F(a 1 (t)) G ft, Vt G 

This means that (/Q, F) : 1Z-^1Z, i = 0,1, and our claim is thus proved. 

Once that the stretching condition in (13. 6p is achieved, the conclusion of the theorem 

follows by Theorem 12.21 11. 

In order to complete the proof, let us verify that any choice of the parameters as in 
(13. 5p and of the domain 1Z = 1Z(xi, yi, z^) in agreement with (Hl)-(Hb) implies that 
conditions (CI), (C2), (C3'), (C4) and (C5) are fulfilled for any path 7 : [0, 1] ->■ K 
joining 1ZJ and 1Z~ lj. In so doing, we will prove that the inequality in (CI) is indeed 
an equality. 

Let us start with the verification of (CI). Since F 3 (x, y, z) = z {\ — ac 3 + ) 

and 7(0) G TZJ = [x e ,x r ] x [ye,y r ] x {z e } = [x e ,x r ] x [y e ,y r ] x {0} by (HI), it then 
follows that 73(0) = and thus = ^3(7(0)) < Zi — 0, as desired. 
In regard to (C2), we have to verify that F 3 \ n - < 0, that is, F 3 (x, y, z r ) < 0, W(x, y) G 
[xi, x r ] xh/i, y r ]. Setting A := x+y, we consider, instead of F^^-, the one- dimensional 
function^ 

4> : [xe + y t , x r + y r ] -)> R, 4>(A) := z r ^1 - ac 3 + 

Computing the first derivative of <\>, we get <p' (A) = z r a(^ J , which van- 
ishes at A = z r . However, since by (H2) we have x^ + y^ > z r , then <fi'(A) < 0, 
VA G [xe + yi,x r + y r ]. Hence, F 3 \ n - < F 3 (xe,ye, z r) and thus, in order to have 
(C2) satisfied, it suffices that F 3 (x£,y£, z r ) < 0. Imposing such condition, we find 

z r (l-ac 3 + ( a i e+ ? l \ 2 ) < 0, which is fulfilled when ^± > . * t+ ™ , 2 . Making z r 




explicit, this holds when z r > y ae "_ 1 (xe + ye) — (xe + ye), that is, when (H2) is 

fulfilled. Notice that the latter is a "true" restriction, since, still by (H2), the right 
hand side of the above inequality is positive. The verification of (C2) is complete. 
As regards (C3 ), we need to check that F 3 (x, y, Ze + Zr ) > % r , W(x, y) G [xe, x r ) x [ye, y r ], 



7 Notice that, by the choice of /Co and Ki, the invariant chaotic set I C fCo U/Ci in Definition 12.21 
lies entirely in the first quadrant and therefore makes economic sense for the application in question. 

8 Just to fix the ideas, in what follows we will assume that 7(0) G TZJ and 7(1) £ 1Z~ . 

9 In several steps of the proof, instead of studying the original problem, through a substitution 
we will be lead to consider a lower dimensional one. Alternatively, we could use the Kuhn- Tucker 
Theorem for constrained maximization problems. We decided to follow the former approach because 
it is more elementary and requires less computations. However, we stress that the two approaches 
require to impose the same conditions (H1)—(H5) on the parameters. 



15 



that is, recalling the definition of S in ( 13.70 . F 3 \s > z r . Notice that, by (HI), 
Ze ~^ Zr = y- Analogously to what done above, instead of F 3 \ s , let us consider the 
one- dimensional function 

z i a A. 
ip : [x t + y e , x r + y r \ ->■ R, <p(A) := 1 - ac 3 + -. 

2 \ ( A + fY 

Since xe + ye > z r > y , by the previous analysis we know that f(A) > (p(x r + y r ) = 
F 3 (x r ,y r , y) . Hence, in order to have F 3 \ s > z r , it suffices that F 3 (x r ,y r , y) > z r , 
that is, 

z r A a(x r + y r ) 

— 1 - QC 3 + —2 | > Z r 

1 V [Xr + yr + f) 

Since z r > 0, making z r explicit, we find 



Zr < 2 f \ -TT( x r + Vr) - ( x r + Vr] 

and this condition is satisfied thanks to (H3). Hence (C3) is verified. 

In order to check (C4), we need to show the two inequalities Fi(x,y,z) < x r , 

W(x, y, z) G TZ and F\(x, y, z) > xe, V(x, y, z) G 71, which are satisfied if 

max Fi(x,y,z)<x r and min Fi(x,y,z) > xe , 

(x,y,z)£lZ {x,y,z)^TL 



respectively 10 !. 

Instead of considering Fi\ n , setting B := y + z and T := [xe, x r ] x [ye + ze, y r + z r \, 
we deal with the bidimensional function 

, - m ™ ^ 2x + B - cAx + B) 2 
$:T^M, $(a;,B):= — -, 

whose partial derivatives are 

d< $> n j 5$ 1 , „, 

— = l-ci(x + B) and — = - - c x (x + B). 
ox oB 2 

Since they do not vanish contemporaneously, there are no critical points in the inte- 
rior of T. We then study $ on the boundary of its domain. 

As concerns $i(B) := §\ {xe}x [ ye+z± y r+Zr ](x, B) = §(x e ,B), we have that $>'i(B) = 
| — C\[xe + B), which vanishes at B = ^- — xe- This is the maximum point of $i if 
B G [ye + ze, y r + z r \. But that is guaranteed by the conditions in (HA). 
Similarly, setting <3> 2 (-B) := ®\{x r }x{y e + Zl ,y r +z r ](x, B) = $(x r , B), we find that its max- 
imum point, still by (HA), is given by B = ^ — x r E [ye + ze, y r + z r ). 



10 



Notice that such maximum and minimum values exist by the Weierstrass Theorem. 
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In regard to $ 3 (x) := ^\[ Xe ,x r ]x{y e +z e }(x, B) = ®(x,yt + ze), we have $ 3 (x) = 1 - 
ci(x + y£ + Zi), which vanishes at x = - — (yi + zA. By the conditions in (HA), x > x r 
and thus $3(2;) is increasing on [a^,a; r ]. Analogously, since x = — (y r + z r ) > x r , 
it holds that $ 4 (x) := §\[ Xe ,x r ]x{y r +z r }(x, B) = §(x,y r + z r ) is increasing on [x e ,x r ]. 
Summarizing, the two candidates for the maximum point of $ on T are (xe, ^ — a^J 
and (x r , — x r ) . A direct computation shows that ^- — a^) < &(x r , ^ — x r ) , 

and thus max iq(x, y, z) = $(x r , ^ — x r ) . Hence, it is now easy to verify that the 

inequality max F\(x,y,z) < x r is satisfied when x r > the latter being among 

the assumptions in (HA). 

The analysis above also suggests that the two candidates for the minimum point of 
$ on T are (xi,ye + ze) and (xi,y r + z r ). Straightforward calculations show that, if 
%l < 5^ (1 - ci(^ + y r + ^ + z r )) , then $(a^, y £ + ^) < $(x e , y r + z r ). Hence, again 
by (HA), min F\(x,y,z) = &(xe, ye + zA. The inequality min F\(x,y,z) > xi 

(x,y,z)£H (x,y,z)£7l 

is thus satisfied when \J Vl ^ 1 — [yt + Z() > xe, which is among the conditions in (HA). 
This concludes the verification of (CA). 

Let us finally turn to (C5). In order to check it, we have to show that 

max F 2 (x,y, z) < y r and min F 2 (x,y,z)>y t . (3.8) 

(x,y,z)£H (x,y,z)£R 

Instead of F 2 \n, setting D := x + z, we deal with the one-dimensional function 

ip : [xe + ze, x r + z r ) — > E, ^(D) := \ D, 

V c 2 

whose derivative is ip'(D) = -^-j== — \. It vanishes at D = which by (H5) is smaller 
than X£ + Z£. Thus max F 2 (x, y, z) = tjj(xe + Z£) and min F 2 (x, y, 2;) = ip(x r + z r ). 

(x,y,z)^Tl (x,j/,z)g7?. 

Hence, the first condition in (13.81) is satisfied if ^(xt+Zf) < y r and the second condition 
is fulfilled if ip(x r + z r ) > y^. It is easy to see that both inequalities are fulfilled thanks 
to (H5) and this concludes the verification of (C5). 
The proof is complete. 



4 Conclusions 



In this paper we have recalled what the SAP method consists in and we have 
applied that topological technique to rigorously prove the existence of chaotic sets 
for the triopoly game model in [13] . By "chaotic sets" we mean invariant domains 
on which the map describing the system under consideration is semiconjugate to the 
Bernoulli shift (implying the features in Remark 12. ip and where periodic points are 
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dense. However, we stress that we did not say anything about the attractivity of 
those chaotic sets. In fact, in general, the SAP method does not allow to draw any 
conclusion in such direction. For instance, when performing numeric simulations for 
the parameter values in ( I3.5p . no attractor appears on the computer screen. The same 
issue emerged with the bidimensional models considered in [13J. The fact that the 
chaotic set is repulsive can be a good signal as regards the overlapping generations 
model therein, for which we studied a backward moving system, since the forward 
moving one was defined only implicitly and it was not possible to invert it. Indeed, 
as argued in [T3], a repulsive chaotic set for the backward moving system possibly 
gets transformed into an attractive one for a related forward moving system through 
Inverse Limit Theory (ILT). In general, however, one just deals with a forward mov- 
ing dynamical system and this kind of argument cannot be employed. For instance, 
both in the duopoly game model in [13] and in the triopoly game model analyzed in 
the present paper, we are able to prove the presence of chaos for the same parameter 
values considered in the literature, except for a bit larger speed of adjustment a. It 
makes economic sense that complex dynamics arise when firms are more reactive, 
but unfortunately for such parameter values no chaotic attractors can be found via 
numerical simulations. 

What we want to stress is that this is not a limit of the SAP method: such issue is 
instead related to the possibility of performing computations by hands. To see what 
is the point, let us consider the well-known case of the logistic map / : [0, 1] — > R, 
f\x) = /ix(l — x), with /i > 0. As observed in [13], if we want to show the presence of 
chaos for it via the SAP method by looking at the first iterate, then we need /i > 4. 
In this case, however, the interval [0, 1] is not mapped into itself and for almost all 
initial points in [0, 1] forward iterates limit to — oo. If we consider instead the second 
iterate, then the SAP method may be applied for values less than 4, for which chaotic 
attractors do exist. Figure 6 shows a possible choice for the compact sets /Co and 
/Ci (denoted in the picture by Iq and I\, since they are intervals) for the stretching 
relation to be satisfied when [i ~ 3.88. This simple example aims to suggest that 
working with higher iterates may allow to reach an agreement between the conditions 
needed to employ the SAP method and those to find chaotic attractors via numerical 
simulations. 

A possible direction of future study can then be the study of economically interesting 
but simple enough models, so that it is possible to deal with higher iterates, in the 
attempt of rigorously proving the presence of chaos via the SAP technique for param- 
eter values for which also computer simulations indicate the same kind of behavior. 
Still in regard to chaotic attractors, we have observed that the SAP method works well 
for models presenting Henon-like attractors, due to the presence of a double folding, 
in turn related to the geometry required to apply our technique. On the other hand, 
a preliminary analysis seems to suggest that the SAP method is not easily applicable 
to models presenting a Neimark-S acker bifurcation leading to chaos. A more detailed 
investigation of such kind of framework will be pursued, as well. 
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Figure 7: The graph of the second iterate of the logistic map with fi ~ 3.88. 



A further possible direction of future study is the analysis of continuous-time economic 
models with our technique, maybe in the context of LTMs, for systems switching be- 
tween two different regimes, such as gross complements and gross substitutes. 

Acknowledgements. Many thanks to Dr. Naimzada, Prof. Pini and Prof. Zanolin 
for useful discussions during the preparation of the paper. 
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